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1 Introduction 



The Lie group method is one of the most powerful method available to analyse 
nonlinear evolution systems and their methods of solution |2 |21 El El- In- 
troduced firstly by Sophus Lie 6 and consists of a systematic procedure for the 
determination of continuous symmetry transformations of a system of nonlinear 
partial differential equations (NLPDE). An important feature of this method is 
that one can derive special solutions associated with nonlinear PDEs straight- 
forwardly which are otherwise inaccessible through other methods. One of the 
obvious advantages of knowing a symmetry group of a system of differential 
equations is that one can use the defining propriety of such a group and con- 
struct new solutions to the system from known ones. One can also use symmetry 
group to explicitly determine types of solution which are themselves invariant 
under some subgroups of the full symmetry group of the system and provides 
a new classification of different classes of solutions. The basic idea of the Lie 
symmetry method is to find the transformations groups associated with a given 
systems under a continuous group of transformations and to find a reduction 
transformation from the symmetries. In particular for nonlinear partial differen- 
tial equations with two independent variables the reduction transformation can 
be used to reduce the number of independent variables by one. In the case of the 
ordinary differential equations the determination of a one-parameter symmetry 
group allows the reduction of the order of the equation by one. Here using the 
criterion of invariance of the equation under the infinitesimal prolonged gen- 
erators, we find the most general Lie point symmetries group of the Thomas 
equation. Looking the adjoint representation of the obtained symmetry group 
on its Lie algebra, we will find the preliminary classification of group-invariant 
solutions. This provides a new exact solutions of the Thomas equation. The 
paper is arranged as follows. In section 2, we introduce the Thomas equation 
and by using the prolongation formulae and the infinitesimal criterion of invari- 
ance we determine the most general symmetry group and the corresponding Lie 
algebra for Thomas equation. Section 3, is devoted to the construction of the of 
group-invariant solutions and its classification which provides in each case new 
exact solutions for the Thomas equation. In section 4, we present conclusions 
and finally in the Appendix, we briefiy give the methods for constructing the 
invariants of the differential equation. 



2 Determination of Symmetry group for Thomas 
equation 

2.1 Thomas equation 

The Thomas equation is an interesting subject in physical sciences [3 HI |n|it 
arises in the study of chemical exchange processus and sometime called Thomp- 
son equation. Proposed by Thomas "7 and has the form: 
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Uxy + + (3Uy + ^UxUy = 0, (1) 

where : = ^,Uy = ^,Uxy = and a, (3 and 7 are constants. For 

exchange processus a > 0, /3 > and 7^6. 

2.2 Lie algebra of Thomas equation 

In order to calculate the Lie point symmetries of the system(l), we first con- 
sider the generator of the group of point transformations. Here there are two 
independent variables x and y, and one dependent variable u. The Thomas 
equation is of second order {n — 2). Then a general vector field on the space 
of independent and dependent variables {X x R where X —< x,y >) takes the 
form : 

V = i{x, y,u)—+ r]{x, y,u)—+ (p(x, V.u) — . (2) 

Where ^, 77 and (/? depend on x, y and u . We wish to determine all possible 
coefficient functions f , 77 and ip, so that the corresponding one-parameter group 
exp ev is a symmetry group of Thomas equation. According the technics using 
in PP, since the Thomas equation is of second order, we need to know the second 
prolongation Pr^v of the infinitesimal generators which is given by the following 
formulae : 

Pr\ ^v + ip-j^ + + (/^""tt^ + '^"^TT^ + (3) 

OUx OUy OUxX OUxy OUyy 

where the coefficients ip^ ^ip^ ^ 'P^^ , ^-^d 9?*'^ are given by : 



= Dxip - UxDxS, ~ UyDxT], (4) 

= Dytp - UxDy^ - UyDyt], (5) 

ip'^y = Dxy{(p - CUx - my) + ^-Uxxy + Wxyy, (6) 

ipVy = Dyyif - 2UxyDyi ~ ^UyyOyT] " UxDyy^^ " UyDyyT], (7) 

(^^^ = DxxV~'^UxxDx^~2uxyDxri~UxDxx^~UyDxxri, (8) 

where Dx and Dy are total derivatives on x and y respectively and Dxy = 
Dx{Dy). After a long calculation, we obtain : 

= 'fix + - £.x)Ux - rjxUy - ^uUl - IJuUxUy (9) 

= (fiy + {ipu - riv)Uy - ^yUx - rjuHy - ^uUxUy (10) 

Vf"^^ = ^xy + iVyu - S.xy)Ux + [Vxu " r]xy)Uy - ^yuU^. (11) 

~l~ {}Puu ^ Cux ^yu^'^x^y ^ux^y ^uu^x^y Vuu^x'^y 

'^^u^x^xy ^^u^y^xy ^y^xx ^x^yy ^u'^y'^xx ^u^x^yy 

+ {y^u - ^x - Vy)'^xy 

^ — ^xx ^" {'^^xu ^xx)'^x 'Hxx'^y '^'^xu'^x'^y (-^^) 



3 



Vu^y'^xx 

— ^yy ~t~ (^'^yu Vyy}^y ^yy^x '^^yu^x^y (■^'^) 
"I" iVuu "^Vyu^Uy ^uuUyUx 

^^u'^y'^Xy ~t~ (^U ^'Uy^'^yy ^T^^^WyWyy f]uu'^y 

^u'^x'^yy 

To use the infinitesimal criterion of invariance, let us introduce A = Uxy + ctUx + 

j3Uy + "fUxUy. 

Suppose that : 

A = 0, if Pr\(A) = 0, (14) 

for every infinitesimal generators v of G, then G is a symmetry group of the 
Thomas equation (1). Substituting the general formulae (9, 10 ) and (11) into 
(3) and (14), after eliminating the dependence between the derivatives coming 
from the equation itself and equating the coefficients of the various monomials 
in the first and second order, we find the defining equations for the symmetry 
group of the Thomas equation to be the following : 



Monomial 


coefficient 


1 


aipx + I^Lfy + Lfxy 


Ux 


Vyu - ixy + IVy - P£,y + arjy 


Uy 


Vxu - Vxy + Ifx - Oirix + Kx 




Vnu - ^ux - Vyu + P^u + ariu + 7Vu 




~^yu ~ "I" Q^^it 




-Vxu - IVx + PVu 


Uiuy 




UxUfi 


Vuu 


^xx 


-^y 


Uyy 


-Vx 


UytlxX 


-in 


UxUyy 


-Vu 



Table 1: Coefficients of various monomials 



The system determining the symmetry group of our equation is : 





afx + f3(py + ipxy 


= 


(15- 


1) 


Vyu 


- ixy + IVy - Piy + aiy 


= 


(15- 


2) 


Vxu 


- Vxy + IVx - oirix + (3^x 


= 


(15- 


3) 


^ux 


- Vyu + P^u + aVu + Ifu 


= 


(15- 


4) 




~£,yu ~ TCj/ "I" Ol£,u 


= 


(15- 


5) 
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-Vxu - IVx + l3Vu = 0; (15-6) 

-^uu = 0; (15-7) 

-Vuu = 0; (15-8) 

-^y = 0; (15-9) 

-Vx = 0; (15-10) 

-^u = 0; (15-11) 

-Vu = 0; (15-12) 



From (15-9) to (15-12) we obtain : 

C = ^(x); and 77 = r]{y) 

Prom (15-4) we deduce that : 

ip{x, y, u) = g{x, y) exp"'^" -\-N{x, y) 

For convenient we can chose "^^^''^^ for g{x,y). 

Using the above results and (15-2) with (15-3) one finds : 

7^^x + P^x = 0; iNy + arjy = 0, 



which yields to 



AT n f ^ 

Nxy = 0; £,x = —, and rjy = —, 



so we have N{x, y) = mx + ny + a where m, n and a are arbitrary constants, 
the Unear first partial differential equation (15-1) leads to am + (3n = 0. 
Then, for ^ — one obtains the most general form of the coefficient 

functions ^, r] and ip which are the form : 



^{x,y,u) = -k-yx + c, (15) 
r]{x,y,u) = k-fy + b, (16) 

V{x,y,u) = -i^^e-^'' + k{(3x - ay) + a, (17) 

7 

where a,b,c and k are arbitrary constants and g{x,y) an arbitrary solution of 
the equation (15-1). Then, we obtain the most general form for vector field v 
which is : 

v = {-k^x + c)^ + {k^y + 6)1- + (-^iM) e-7« + k{px -ay) + a)^. (18) 
ox ay 7 au 
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Hence, the Lie algebra of infinitesimal symmetry of the Thomas equation is 
spanned by the four vectors fields : 



V3 



d_ 

d_ 
dy' 

d_ 
du 



and the infinite dimensional subalgebra : 



9{x,y) ^ 
7 



-7U 



du 



(19) 
(20) 
(21) 

(22) 
(23) 



The one-parameter groups Gi generated by the Vi are given in the following 
table. The entries gives the transformed point exp{evi){x,y ,u) = (x,y,u) : 



Gi : {x + e,y,u) 
G2 : {x,y + e,u) 
G3 ■■ {x,y,u + e) 

G4 : (xe-^%ye''%^a;(l-e-^^) + -2/(l-e^^)+w) 
7 7 

Gg : {x,y,-\og['yg{x,y)€ + e''^]) 
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(24) 
(25) 
(26) 

(27) 
(28) 



Since each group Gi is a symmetry group, then if u = /(x, y) is a solution of 
the Thomas equation, so are the functions : 

e,y) (29) 

-e,) (30) 

) + e (31) 



Ul 


= /( 


U2 


= /( 


U3 


= /( 






U4 










1 


Ug 
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a 
—1 

1" 



(32) 
(33) 



where e is any number real. 

The commutation relations between these vector fields are given in the following 
table : 
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[Vi,Vj] 




V2 


V3 


V4 




Vl 











-jvi + /3v3 




V2 











7^2 - av3 




V3 














'"-If 


Vi 


7Dl - I3V3 


-7D2 + aua 














-Vfy 










Table 2: Commutator table for the Lie algebra vi and Vg 



where 

ij} = -7a; + '^yfy - 7(/3a; - ay)f. 

Since the totaUty of infinitesimal symmetries must be a Lie algebra, we can 
conclude that if f{x,y) is any solution of the equation in (15-1), so are f^, fy 
and --fxf^ + jyfy - ^{f3x - ay) J. 

3 Classification of Group-Invariant Solutions 

Recall first that in general to each one parameter subgroups of the full sym- 
metry group of a system there will corresponding a family of solutions, such 
solutions are called invariant solutions 1 . In this paper we are interesting only 
on the symmetry algebra g of the Thomas equation which is spanned by the vec- 
tor fields Vl, W2, V3 and Vj^. Every one-dimcntional subalgebra of g is determined 
by a nonzero vector v of the form : 

V — aivi + a2V2 + a^vz + OiVi, (34) 

where ai are arbitrary constants. Our task is to simplify as many of the coeffi- 
cients ai as possible through application of adjoint to v. To compute the adjoint 
representation, we use the Lie series : 

Ad{exp{evi))vj = vj - e[vi, vj] + y [u^, [vi,Vj]] (35) 

Then from the commutation table 2, we obtain the following table : 



Ad 


Vl 


V2 


V3 


V4 


Vl 


Vl 


V2 


V3 


V4 + e7Wi — ef3v3 


V2 


Vl 


V2 


V3 


V4 — t'^V2 + eav3 


V3 


Vl 


V2 


V3 


V4 


Vi 


e-i^vi + ^il-e-^')v3 


e^'v2 - ^{e-^' - 1)V3 


V3 


V4 



Table 3: Adjoint table for the Lie algebra Vi 
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In what follows, we begin the classification process 
Case.l 



Suppose first that 04 ^ 0, we can assume that a4 = 1, and the vector field v 
takes the form : 

V ~ oiwi + 02^2 + 03^3 + ^^4 (36) 

Refereing to the table. 3, if we act on such a by Ad(ex'p[^)v) , then we can 
make the coefficient vanish . Then the reduced vector field takes the form : 

d d d 

v = {ax- ix)— + (a2 + 72/)^ + {Px - ay)^. (37) 

The invariants <; and x can be found by integrating the corresponding charac- 
teristic system (see Appendix for construction of the invariants), which is : 



dx dy du 



ai — 7a; a2+ 72/ Px — ay 
The obtained solution are given by 



(38) 



X = (ai - ix){a2 + 72/), (39) 

and 

, = u+^-x+ i^^l±^/o,(7x - a,) - , (40) 

7 7^ 7^(70; — ai) 

Therefore, a solution of our equation in this case is : 

u = f{x, X, '?) = '? - { ^^^ \"'°'^ )log{'yx - ai) + „, — - -x (41) 

Y Y[lx — aiY 7 

The derivatives of u are given in terms of c and x as : 

= -7(o2+72/)Tx 7 7 ; (42) 

7(7X - ai) 7 

= -- +7(ai -7a;)?x; (43) 

7 

Wxy = -7^'?x - 7^<?xx; (44) 

Substituting (42), (43) and (44) into the Thomas equation we obtain reduced 
equation of the form : 

7^XCxx + + 7^(7 - - aa2)Cx + y = 0- (45) 

For Cx = ^, one obtains the following equation : 

e' = -,e^ - (IZ^^I^^ + ^ = 0; (46) 
7 X 7 X 
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this Riccati type equation is transformed by using the change : 

y{x) = exp( j 7^(X)), (47) 
multiplying by x one obtain the Fuchs type equation : 

X^y" + eX2/' + fnxv = 0, (48) 
where e = ^^"^"^"""^^ and m = ^. 

This equation admits a solution developable in series in the neighboring of 
zero, for X > 0, we have : 

oo 

2/p(x) = E'^»^"' (^9) 

n=0 

where ao = 1 and a„ = , n > 1 

In the case x < by introducing the transformation x i — > —X) This series 
have an infinite radius of convergence for all x > 0) then if we introduce the 

change Zp{x} = Zp is a particular solution of the equation (46). Putting 
9 = z + Zp, (46) becomes a Bernoulli type equation : 

z = —70 — 2^ZpZ, (50) 
which is transformed in the linear equation (with z = j) 

/' =27^ + 7. (51) 

The solution of (51) is : 

fix) = 9p{x)yl (52) 

where Qp ia a primitive of ^ 
Therefore, 

^^(x) = -^+^p(x)- (53) 

Thus, we have c = ?(x) = / 0{x)dx, 

and then, we obtain the solution in this case which is : 

u = <;ix) - { ^^^ ^'^^^ )log{'yx - ai) + ^ - (54) 

Y 7^(70; — ai) 7 

where X = («i - 7a;) (^2 + 72/)- 
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Case.2 

Suppose now that 04 = 0, and 03 ^ 0. In this situation the invariants asso- 
ciated to the field vector v are of the form : 



Case 2.1 : If 0102 ^ 0, 



V = aiwi + 02^2 + W3; (55) 
d d d 

ox ay ou 



The invariants are : 



X = a2X-aiy\ (56) 

y_. 

02' 



? = ; (57) 



Next, we have : 



Ux = a2Cx; (58) 

1 , , 

Uy = aiCx; (59) 

02 

Wxy = -aia2Cxx- (60) 

Substituting these quantities in the Thomas equation, we find the reduced equa- 
tion : 

- aia2Cxx + {010.2 - (3ai + 7)^^ - aia27^Y H =0; (61) 

a2 

Putting = 0, then (61) takes the form : 

- 01026*' + (aa2 - Pai + 7)6* - 01027^^ + — = 0; (62) 

(12 

which is the Riccati type equation. 

Case 2.1. a : If {aai - pai + jf + Aai^(3 > 0. 

Then, our equation admit one constant solution 60 which satisfying the sec- 
ond order equation: 

(aai — (3ai + — 010272:^ + — = 0. (63) 

02 

Then by redefining the variable = f + 9o, one obtain the following Bernoulli 
type equation : 

- aia2f' + {aa2 - /3ai + 7 - 201027^0)/ - aia27/^ = 0; (64) 
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Put now f = \ one obtain a linear equation: 

, 2aia2'yOo -aa2+ pai+j 

h = h + -f (65) 

aia2 

This equation admit a solution, since ■^»i»2'yeo-aa2+0a,+j ^ ^^^^-^ . 

2aia27^o - aa2 - + 7 ^ 7^1 ^2 

h = Aexp{ x) - 2 3 — — (66) 

aia2 2ai 0276*0 - aai — pai + 7 

where A is an arbitrary constant. Hence, 
C = / o a J' dx + Oox + cte (67) 

J ^ex/y^^ ^^^^ ^; 2aia2j0o-aai-l3ai+"/ 

After integration we obtain : 

<;=-log{A-^e-''^) + OoX + cte (68) 
7 G 

where C = 2aia27go--a^a2+/3ai-7 _ Finally, the solution of the Thomas equation 
is: 

U = -logiA - le-Cia2X-ary) ^ g^^^^^ _ ^^y^ _^ y_ ^ (gg) 

70 a2 



Case 2.1.h : If {aa2 - Pai + 7)^ + 4ai7/3 < 0. 

In this situation wc have : 

(aa2-/3ai+7) ^_ 2^^^0_ ^^0) 
aia2 aia2 

The reduced equation obtained can be written in the separation form as : 

de 



{aa2-0ai+l) _ Q2 i 
aia2 ^ aia'i 



dx. (71) 



If wc put Ai = /3ai+7) , _ _ _ — thus the problem becomes 

to the integration of the following equation : 

Let S = ^^^^T^, S > 0. Then, we look to the integration of 



1 
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By using the change of the variable t = 'E+ , we have : 



2A; 

and therefore, 



f d6 1 .^+^^ 



then, 



and we obtain : 



1 

—z=arctang{ 7=^) = -^2X + cte, (75) 



= VEtg{A2VEx + cte) - (76) 

2A2 



= VEtgiA2VEx + cte)-^^, (77) 
= VeJ tgiA^VEx + cte)dx - + cte, (78) 

= ^^05(l + t5((^2\/Sx + cte)- ^X + cte. 
Next, we find the solution of the Thomas equation which is : 



u = -^log{\ + tg{Ai\pB.{a-2.x - aiy) + Aq) - 7^(022; - aiy) + — + cte; (79) 
where Aq = constant. 
Case 2.2 : If 02 = 0, ai 7^ and oi 7^ ^ 
In this case the vectors field is: 

^ = «iijr + ir:- (80) 



The invariants are : 



dx du 



X = y, (81) 

<r = ~aiu + x, (82) 



then, u = ^ + and 

' ai ai ' 



(83) 
(84) 

0. (85) 



ai 

Uy = ^, (84) 
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After substituting these expressions into the Thomas equation, we obtain the 
following reduced equation : 



and then, 
Next we have, 
and therefore, 



(/3oi +7)<rx = aia; (86) 



^ = -Z^^X + cte. (87) 



/3ai +7' 



^ " C + cie; (88) 



ai /3ai + 7 

X a 
ai (3ai + 7 



u = — - — y + cte. (89) 



Case 2.3 : If ai = ^. 

The vector field v take the form : 



and the invariants are 



Therefore 



and then. 



—7 d d 



X = y; (91) 

<j = /3a; + 7U. (92) 



u = ^ - -X, (93) 

7 7 



7 
7 

0. (96) 



(94) 
(95) 



Substituting the above expressions into Thomas equation, one find that a/3 = 0. 
Then we conclude that there is not the solution in this case. 

Case 2.4 ■ If ai = 02 = 0. 

we have : 

v = l: (97) 
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the invariants arc <; = x, and X = J/, then hke in the last case there is not the 
solution of the equation. 



Case.3 



In this case we choose 04 = 03 = 0, then the vector field is the form : 



Case 3.1 : 0102 7^ 0. 

The vector field is of the form : 



9 d 



and the invariants <j et x are solutions of the following characteristics system : 

(100) 



dx dy du 



The solutions are 



X = x-^; (101) 

(12 

<; = u. (102) 
next we have : u = <r, therefore : 

Ux = Cx; (103) 

Uy = (104) 



a2 

02 



(105) 



By substituting these expressions into the Thomas equation, we obtain the 
following equation : 

putting 6 = <;^, we obtain the Bernoulli equation : 

9' + a2a)0 + 76*2 = (107) 
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Case 3.1. a : If 02 = §. 

The reduced equation takes the form : 

e' + 76*2 = 0; 

which have the solution : 

0= ' 



IX + cte 
and therefore, 

1 y 

u = —log{'-f{x ) + ko) + cte; 

7 0,2 

where fco, is arbitrary constant. 
Case 3.1.b : If 02 9^ ^. 

We adopt the change = ^, we obtain the following linear equation 

z' = {P- a2a)z + 7, 
which admit the general solution : 



7 



/3 — 02a 

Therefore, 



p — a2a 

and then, 

-1 7 

u = — ^05(1 H 7^ r-, — ^-7 ) + cte; 

7 k{p — a2a)e^ 02'— 7 

where A; = is arbitrary constant. 

Case 3.2 : If ai = 0. 

In this situation we takes v as : 

d 

oy 

and the invariants are : 

X = X 
c; = u 
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Then u — and therefore, 



Uy = 0, 

Uxy = 0. 



(118) 
(119) 
(120) 



Substituting these expressions into the Thomas equation, then the sokition of 
the corresponding equation is a constant : 



4 Conclusion 

In this paper by using the criterion of invariance of the equation under the 
infinitesimal prolonged infinitesimal generators, we find the most Lie point sym- 
metries group of the Thomas equation. Looking the adjoint representation of 
the obtained symmetry group on its Lie algebra, we have find the preliminary 
classification of group-invariants solutions. We have seen that the obtained re- 
duced equation in such case can be transformed on known equation by using 
an appropriate change of the variables. It is intersecting to extended this con- 
struction to the supersymmetric Thomas equation, more detail in this question 
will be given elsewhere ^13,. 




(121) 
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5 Appendix 

In this Appendix, we introduce the method for finding the invariants of the 
group generated by the foUowing infinitesimal generator : 

d d d 

v = {ai- jx)— + (02+72/)^ + {I3x - a^)^- (122) 

The invariants are found by integrating the corresponding characteristic system, 
which is 

dx dy du 



ai — 7X 02 + 72/ l3x — ay' 
The first of these two equations : 

dx dy 
ai -jx a2 + jy' 
is easily solved, the solution are : 



(123) 



(124) 



- log(ai - 7a;) = - log(a2 + 72/) + c. (125) 
7 7 

So, one of the invariant is : 

X = {ai-'yx){a2 + 'yy). (126) 

Note that x is a constant for all solution of the characteristic system, so we can 
replace y by — ^ before integrating. This leads to the equation : 



dx (I I 



(127) 



ai - 7X Px- , "X , + ^ ' 

^ ' ^ 7(ai— 7x) 7 

therefore : 

{— 7 + ^--}dx = du. 128 

ai — jx 7(ai — •yxy 7(01 — jx) 

Then the equation has the following solution : 

- (^^1±^ iog(a, - 7.) - =u + k, (129) 

7 7^ 7 (di — 7a;) 

for k an arbitrary constant. Finally we have : 

/3 {/3ai + aa2) ax 

? = w+-a;H 2 log(ai - 7a;) + -, (130) 

7 7^ 7 (ai — 7a;) 

is the second invariant. 



17 



References 

[1] Olvcr P.J, Application of Lie groups to Differential Equation, Springer, 
New York 1986. 

[2] Stcphani H, Differential equations: Their Solution Using Symmetries, Cam- 
bridge University Press, 1989. 

[3] Bluman G.W and Kumei, Symmetries and Differentials Equations, 
Springer, New York 1989. 

[4] Ibragimov N.H, CRC Handbook of Lie Group Analysis of Differential Equa- 
tion, CRC Press, Boca Raton, 1986. 

[5] Rodriguez M.A. and Winternitz. P, J.PhysA:Maths.Gen 37(2004) and the 
references quoted therien. 

[6] Lie S and Engel F, Theories der Tranformation gruppen Vol.39 (Leipzig 
Teubner) 1890. 

[7] Wong W.T. and Fung P.C.W, Linuo cimento Vol.99B 1987.163. 

[8] Zheng K et al., Physica Scripta Vol.40 1989, 705. 

[9] Sakovich S.Yu, Journal Physiacs A: Maths.and Gen.21 1988,L1123 
[10] Guang-Mei.W et al., Gechoslovok Journal of Physics Vol.52 2002, 749. 
[11] Zhenya Y, Gechoslovok Journal of Physics Vol.53 2003, 298. 
[12] Thomas H.C Journal American Chemistry Soc.66 1944, 16 64. 
[13] A.Ouhadan and E.H . El Kinani, in preparation. 



18 



